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Abstract. We study implications of unitarity for pseudo-orbit expansions of the 
spectral determinants of quantum maps and quantum graphs. In particular, we 
advocate to group pseudo-orbits into sub-determinants. We show explicitly that the 
cancellation of long orbits is elegantly described on this level and that unitarity can be 
built in using a simple sub-determinant identity which has a non-trivial interpretation 
in terms of pseudo-orbits. We reformulate Newton identities and the spectral density 
in terms of sub-determinant expansions and point out the implications of the sub- 
determinant identity for these expressions. We analyse furthermore the effect of the 
identity on spectral correlation functions such as the auto-correlation and parametric 
cross correlation functions of the spectral determinant and the spectral form factor. 
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1. Introduction 
1.1. Overview 

When calculating quantum spectra with the help of periodic orbit sums such as, for 
example, arising from semiclassical expressions, one typically encounters problems due to 
divergencies resulting from summing over a large number of periodic orbits which grows 
exponentially with length. This is in particular the case for quantum systems whose 
underlying classical dynamics is chaotic pQ. To apply these periodic orbit expressions 
for determining quantum spectra, the number of relevant orbits needs to be reduced. 
This is either achieved by reordering the orbit contributions making use of cancellations 
such as is done in the cycle expansion [2] or one can also utilise unitarity of the quantum 
dynamics leading to additional relations between the coefficients of the characteristic 
polynomial and thus to finite sums over pseudo-orbits [3], SI El EJ, [Zl E] ■ 

A related problem is the semiclassical calculation of spectral correlation functions. 
They are conjectured to follow Random Matrix Theory (RMT) for quantum systems 
with chaotic classical limit. Establishing this connection explicitly using semiclassical 
periodic orbit formulae for the spectral form factor could only be achieved fairly recently 
following the work in [9]. This calculation has been extended in [10] yielding the full 
spectral form factor as predicted by RMT for times smaller than the Heisenberg time 
Th- (This is the time needed to resolve distances of the order of the mean level spacing 
in the Fourier-transformed spectrum) . The spectral form factor for times larger than T# 
has been obtained using semiclassical periodic orbit expressions in [TT]. The calculation 
is based on a generating function approach containing two spectral determinants both 
in the numerator and denominator at four different energies. The derivation makes ex- 
plicit use of the fact that the spectral determinant is real for real energies. Although this 
is obvious from its definition, Eq. (|8j below, it is not clear a-priory when considering 
the representation of the spectral determinant containing periodic-orbit sums. A real 
spectral determinant in terms of periodic orbits can only be semiclassically obtained by 
exploiting periodic-orbit correlations due to unitarity. 

The above problem illustrates that we need a better understanding of correlations 
between periodic orbits and in particular correlations between long and short orbits. 
To analyse these correlations in more detail, we study here quantum unitary dynamics 
described in terms of finite dimensional unitary matrices. We will in particular advocate 
to consider spectral quantities in terms of sub-determinant expansions. We will derive 
sub-determinant expressions for a range of important spectral quantities and consider 
these for examples such as quantum graphs [12]. 

The paper is structured as follows: We first introduce the spectral determinant and 
explain the known implications of unitarity for this quantity. We analyse in Sec.|2]further 
implications of unitarity on pseudo-orbit expansions. In this context, we present a sub- 
determinant identity for unitary matrices and we discuss the implications of this identity 
on Newton identities and a pseudo-orbit expansion of the spectral density. In Sec. |3j we 
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derive expressions for spectral correlation functions such as the auto-correlation and the 
parametric cross correlation function of the spectral determinant and the spectral form 
factor in terms of sub-determinant expansions. Implications due to the sub-determinant 
identity will be discussed. 

1.2. Some basic properties of the characteristic polynomial of a matrix 

Consider a general complex matrix U of dimension N . Its characteristic polynomial is 
given by 

N N 

Pu{z) EE det (Z-U) = J2(-i) N ~ n aN-nZ n = HO - Z nl i 1 ) 

n=0 n=l 

where the complex numbers z n are the eigenvalues of U. The complex coefficients a n of 
the polynomial in (JTJ will be at the centre of interest in this article. Here, do = 1 and the 
remaining N coefficients a n , n = 1, . . . , N, are N complex numbers which contain the 
same information as the N eigenvalues z n . Note that the characteristic polynomial is 
invariant under conjugation U h- > CUC^ 1 with a non-singular matrix C. The coefficients 
a n are thus matrix invariants (as are the eigenvalues) and can be expressed in terms 
of other matrix invariants such as traces of powers of U. Indeed, expressions for the 
coefficients a n in Eq. (JTJ) in terms of eigenvalues or traces of U can be easily written 
down, for instance 



N 



^2 z n = tr U, 



n=l 

a-2 = g ZnZm = 2 ^ tr2 U _ tr U ^ ' 

Similar formulae expressing the a n 's in terms of traces hold for all n [13J. Note however, 
that a at = Yln=i z n = det U has a much simpler expression in terms of the determinant 
of U. 

Alternatively one may express the coefficients in terms of sub-determinants of U. 
Denote the set X = {1,2, .. . ,N} and let T C X be some subset of X of cardinality 
|T| < N. Then T defines a quadratic |T| x |T| submatrix Ur which is obtained from U 
by keeping only those rows and columns with indices belonging to T. We will denote 
the determinant of Ur as 

d r = det U r . (2) 

Using the linearity properties of the determinant with respect to its rows (or columns), 
it is then straight forward to show that 

a n = dr . (3) 

TcX: |T|=n 

The sum extends over the ( ) different choices of n rows (and the corresponding 
columns) that build the sub-matrix Ur- While a n is a matrix invariant it is noteworthy 
that this is in general not the case for the individual contributions det Ur- 
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2. On pseudo-orbit expansions in terms of determinants 

2.1. Basic relations 

Let us now consider the characteristic polynomial and some related spectral functions 
for the specific case of unitary matrices U. We will keep the discussion general here and 
will only later refer to U as the evolution matrix for a quantum system. 

A unitary matrix U of dimension N has N uni-modular eigenvalues z n = e t9n . This 
implies the functional equation 

Pu(z) = (-z) N e^Pu(l/z*y (4) 

for its characteristic polynomial where z* denotes the complex conjugate of z and 
e %4> = ^g^- jj — aN _ Comparing coefficients of z n on both sides of the functional equation 
Q results in the explicit relation 

a N _ n = e*a* n (5) 



between the coefficients of the characteristic polynomial. Below in Sec. |2.2[ we will 
generalise this relation to individual determinants of sub-matrices contained in the 
coefficient a n according to (j3j). 

For unitary maps it is useful to introduce the following variant of the characteristic 
polynomial, the so called zeta function, 

N 

Cu(0) = e-* m Pu{e w ) = det (l - e~ i6 U) = ^(-l)"a„ e -^ . (6) 

n=0 

This is a 2ir periodic function in the variable 9 which vanishes exactly at the spectrum 
of real eigenphases {9 n }^ =l . In terms of sub-determinants ^ one may also write 

Cu(9) = J2dre^ 9+ ^ (7) 
r 

where the sum is over all subsets T C X including the empty set T = with |T| = |0| = 
for which we set d$ — 1. 

The functional equation Q implies that 

Z V {6) = e^-*ku(9), (8) 

usually referred to as spectral determinant, is real for real 9, i.e. Zjj{9)* = Zu(9). 
Another spectral function which will be discussed later is the density of states 

N 

p{9) = Y,^{9-9 n ) (9) 

71=1 

where 5 27r (x) = Y^=-oo $( x + 27rn) is the 2-n periodic 5-comb. The density of states can 
be expressed as 

p(9) = -^-lm\ogZ u (9 - ie ) = ^- + IAi m i ogC[/ (0 - ie ) (10) 
7r au 271 71 da 

in the limit e — > 0. This expression directly leads to the trace formula which expresses 



the density of states in terms of periodic orbits. We will discuss this in Sec. 2.3 together 



with a novel expansion in terms of sub-determinants presented in the next section. 
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2.2. A sub -determinant identity for unitary matrices 

We will now give an identity between the sub-determinants dr for unitary matrices which 
contains much more detailed information than Eq. ([5]). We will interpret this identity 
in terms of periodic orbits and will discuss the implications for spectral measures in the 
remainder of the article. 

Theorem: Let U be a unitary matrix of dimension N with determinant det U = e %< ^ and 
T C X = {1,2, N} with n = \T\. Denote the complement of T in X by f = X — T . 
Then the following identity for the determinants of the submatrices Ur and Up holds: 

det U T = (det Up)* . (11) 

Proof. The identity can be proven directly by writing U, without loss of generality, 
in block-form 



U= \ Z rr ( 12 ) 




where V is an n x (TV — n), and W an (N — n) X n matrix. Let us first consider 
the case where the n x n block Ur and the (N — n) x (iV — n) block Up are both not 
singular. One then has e 1 ^ = det U = det Ur det (U T — WU^V) . Unitarity of U implies 
U^VUt = -Wt an d U r Ut + WW^ = l {N _ n)x{N _ n) such that 



which proves the identity (11) for det Ur ^ and det Up ^ 0. Continuity of det Ur 
allows us to extend the identity to singular blocks Ur and Up. As a corollary Ur is 
singular if and only if Up is singular. ■ 

Given that this is a very basic property of unitary matrices and that the proof is 
fairly straight forward, it would surprise us if this identity had not been noted previously, 
although we do not know of an explicit reference. It implies some fundamental 
connections between orbits and pseudo-orbits of dynamical systems, which, in our view, 
are worth exploring. We will discuss these implications in the following sections. 

As a straightforward consequence, one obtains for the zeta function ^ for iV odd, 

Cc(0)= J2 {dre~ l ^ m + d r e^e-^ 9+ ^ N -^) . (14) 

T:|r|<V/2 

The formula remains true for iV even if appropriate care is taken for contributions with 
|T| = N/2; only half of these contributions should be counted and this half needs to 



be chosen appropriately. Expression (14) resembles Riemann-Siegel look-alike formulae, 
see [3111]. 

2.3. Pseudo-orbit expansions in terms of determinants 

In the previous sections, we have expressed the characteristic polynomial Pu(z) and 
related expressions in terms of the determinants dr- Before we turn to express the 
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density of states or spectral correlation functions in a similar fashion, we will consider 



how the identity (11) can be interpreted in a periodic orbit language. To this end, 
we briefly explain what we mean by a 'periodic orbit' in terms of a finite matrix and 
introduce some related notation. Analogous finite pseudo-orbit expansions in terms of 
short orbits have recently been discussed in the context of quantum graphs [8]. We 



stress here expansions in terms of sub-determinants which together with Eq. (11) give 
compact expressions for spectral quantities in terms of short periodic orbits. 

2.3.1. Periodic orbit representations In the present setting of a unitary N x N matrix 
a periodic orbit p = i\ . . . i n of (topological) length \p\ = n is a sequence of n integers 
i m G {1, 2, . . . , iV} where cyclic permutations are identified, e.g. 134 = 341. One should 
think of a periodic orbit as a set of indices of the matrix U that are visited in a periodic 
way. A primitive periodic orbit is a sequence p — %\ . . . i n which is not a repetition of 
a shorter sequence. If p is not primitive we denote its repetitions number by r p . An 
irreducible periodic orbit never returns to the same index, that is, all i m are different; 
the length of an irreducible orbit is at most N. We also define the (quantum) amplitude 

n 
m=l 

of a periodic orbit. If p is not irreducible one may write its amplitude as a product of 
amplitudes of irreducible orbits, for instance tj2l3 = %2%3- 

A pseudo-orbit 7 = p^p™ 2 ■ ■ ■ p™ n with non-negative integers m ; is a formal product 
of periodic orbits pi with length I7I = J2i m i\Pi\ an d amplitude t 7 = O- We will say 
that a pseudo-orbit is completely reduced if it is a formal product of irreducible orbits 
and irreducible if all mi are either one or zero and if any given index appears at most in 
one pi with m\ — 1. 

These definitions allow us to write the trace tr U n = Y^p-\ p \= n F^p as a sum over 
amplitudes of periodic orbits of length n. Using logdet(l — e~ l6 U) = trlog(l — e~ l9 U) 



in (10) and expanding the logarithm one arrives at the trace formula 

peP r=l 

where sum over p extends over the set of all primitive orbits denoted by V, and the 



additional sum is over all repetitions. Here, like in Eq. (10), it is always understood 
that 9 = 9 — ie and the limit e — > is taken. 

Performing the sum over repetitions in the trace formula shows that it is equivalent 
to an Euler-product type expansion 

Cu(9) = II (! - tve- iW ) ■ (17) 

Note that this is an infinite product which only converges for sufficiently large e > 
and analytical continuation is necessary to move back to the axis e = 0. Such an 
analytic continuation is of course given by the expression (J6]) which is by definition a 
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finite polynomial in z = e l9 . Strong correlations between the amplitudes of long and 
short periodic orbits have to exist in order to reconcile both expressions. Indeed, large 



cancellations can be shown to exist by expanding the product (17) and ordering the 
terms with increasing orbit length such as the cycle expansion proposed in [2]. After 
expressing amplitudes of reducible (arbitrarily long) orbits as product of amplitudes of 
irreducible (and thus short) orbits, the cancellation mechanism emerges [6l[7j. 

Revisiting Eq. ([7]) and observing that each determinant dr can indeed be written 
as a sum of |T|! irreducible pseudo-orbits 7 of length I7I = |T|, we obtain 

^r=£(-ir +1 V (18) 

7 eP r 

Here, Pr is the set of all irreducible pseudo-orbits which cover the set T completely, 
that is, which visit each index in T exactly once. There is a one-to-one correspondence 
between these irreducible pseudo-orbits and permutations. Indeed any permutation of 
symbols in T can be written uniquely as a product of cycles such that each symbol 
appears exactly once (up to the ordering of the cycles which is irrelevant as they 
commute). Each such product of cycles, that is, each irreducible pseudo-orbit, defines a 
unique permutation. We denote the number of cycles (irreducible orbits) that make up 
a given pseudo-orbit 7 as 07 such that (— 1)°"t +1 gives the parity of the permutation. 



2.3.2. Interpretation of the identity (11) in terms of periodic orbits Everything said in 
the previous subsection is valid for general, not necessarily unitary matrices. Unitarity 
leads to further non-trivial relations between the amplitudes of short and long orbits 
such as the functional equation Q resulting in the relation ^ for the coefficients of the 
characteristic polynomial which can in turn be written in terms of orbits. 



In Sec. 2.2 we showed that there is a much more detailed link between sub- 
determinants and thus orbits. The identity (11), dr = e^d^,, also provides a connection 



between short and long orbits, but it has in addition an interesting interpretation 
in terms of linking irreducible pseudo-orbits in different parts of 'phase space'. T 
and its complement T are by definition disjoint and its union forms the whole set 



X = {1, 2, . . . , N}. As stated in Eq. (18), dr,d*^ consist of all irreducible orbits and 
pseudo-orbits which completely cover the set T, T, respectively, (passing through every 



index in each of the sets exactly once). The relation (11) thus implies that the sum 
over all irreducible pseudo-orbits that cover V is equivalent in weight to the sum over 
all irreducible pseudo-orbits that cover its compliment T. The two contributions from 
the pseudo-orbits in T and the complement T yield together a real term in the spectral 
determinant, as the contributions from T and T are complex conjugated to each other up 
to a global phase. This has direct implications for periodic-orbit approaches to spectral 
statistics in quantum graphs and quantum maps in the semiclassical regime, as discussed 
in Sec. 13.31 

Given the close relationship between unitary matrices on the one hand and quantum 
maps and continuous quantum systems on the other hand, this finding has potentially 
far reaching consequences. It suggests that the quantum (and semiclassical) weights 
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associated with periodic orbits and pseudo-orbits of classical maps and flows are spatially 
correlated on all levels. In particular, summing over all orbits associated with a given 
subset of the full phase space should yield a total amplitude which is equal to the 
contribution from the orbits in the complement and both contributions are phase related. 
We will leave it at these general remarks here and restrict the discussion again to finite 
dimensional unitary matrices only. 



2.3.3. Density of states and Newton identities We will now consider the density of 
states and show that it can be expressed in terms of completely reduced pseudo-orbits. 
Equivalently, one can write it as a sum over products of sub-determinants dr- The 
latter has the advantage that these expressions keep track of the relation (11) between 
individual determinants which is lost on the level of pseudo-orbit sums. 

Making use of Eq. ( 10 ) , we would like to express log (u in terms of sub-determinants. 
We do this by exploiting the identity 

i 00 r-2-K 

- log(l - x) = — J> - 1)! / e Mxe ~ ta da , (19) 

71 n=l ^° 

which formally requires x < 1. Setting 1 — x = Y2r dre~^ e+7T ^ and using ([7]), we 
formally obtain \og(u on the left hand side of Eq. ( Jl9| . After expanding out the 
exponentials, interchanging integration and summations and carry out the integration 
over a, one obtains 

\og(u(9) = - (H - l)!e- i(9+ " )|mrH " |m| Y[ . (20) 

m:|m|>0 j 3 

Here m = (mi, . . . , m 2 N_ 1 ) is a set of 2 N — 1 non-negative integers and Fx, . . . , T 2 n_i is 
some enumeration of all non-empty subsets T C X. The integer rnj is the multiplicities 
of that subset Tj in one contribution to (20). We have also introduced the notations 
|m| = Yl^i 1 m j an d |mr| = Y^j=\ l m j\Fj\- Note that an analogous equation to (20) 
can be given either in a coarser way in terms of coefficients a n or in a more detailed 
way in terms of products of irreducible pseudo-orbits. The expression in terms of the 
determinants dr is the most detailed one in which the relation (11) between long and 
short orbits remains explicit. 

Before moving on to the density of states let us consider the well known expansion 
— log £[/(#) = J^^Li \e~ t6n tr U n and compare the coefficients of e~~ tdn with the 
corresponding ones in (20). This gives us a direct way to express the n-th trace in 
terms of the sub-determinants <ip, that is, 

trl7» = n(-l) B d m l " X ) ! II ^f" • ( 21 ) 

m:|mr|=n j=l J 

This formula is reminiscent of the well-known Newton identities that express the traces 
of powers of a square matrix in terms of the coefficients of the characteristic polynomial, 
see for example [14] . Indeed, as mentioned above, there is an expression of the form (20) 
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in terms of the coefficients a n instead of the dr- The corresponding derivation of the 



traces leads to the Newton identities. In (21), we have in fact derived a more detailed 



identity; it allows us to express the (arbitrarily long) periodic orbits that add up to the 
traces tr U n explicitly in terms of pseudo-orbits of length smaller than the matrix size 



N. Furthermore, using (11), one has an explicit expression of traces of any power in 
terms of pseudo-orbits of maximal length N/2. Ordering the sequence (Ti, . . . ,r 2 iv_i) 
such that it is non-decreasing in length, then 



>V " 1 -- 1 (-d r .) m i Cr 1 (-^.e^n- 



trl7» = n(-l)» £ (H-l)! J] ~^th II m , ( 22 ) 

m:|mr|=n j=l 3 i=2 JV - 1 3 

gives the n-th trace in terms of contributions which can be computed from irreducible 
orbits of length smaller than N/2. In the following it will always be understood that 



products of the form appearing in (21) may be expressed analogously to (22) in terms 
of short orbits. 



Eventually the density of states follows directly from (20): 



p (0) = ± - Im J2 (|m| - iy)j^ e -^)\™r\-™\™\ J] _£z_ . ( 23 ) 

m:|m|>0 j 3 

3. Spectral fluctuations in terms of sub-determinants and short orbits 

There is a wide variety of measures for spectral fluctuations which have been considered 
in the past. We will focus here on expressing spectral measures in terms of 



sub-determinants and show how the relation (11) can be used to understand the 
contributions of long orbits. We will in particular consider ensembles of unitary matrices 
where the ensemble average corresponds to an average over system parameters or 



disorder. In Sec. 3.3 we will also discuss applications which only involve a spectral 



average for a fixed physical system. 
3. 1 . Spectral fluctuations 

For a given ensemble of unitary matrices we denote the ensemble average of some 
quantity f(U) by (f(U))u- In the following, we will consider cross- or auto-correlation 
functions for the spectral determinant, the density of states and other quantities. We 
start by giving some general definitions. 

3.1.1. The auto- correlation function of the spectral determinant This auto-correlation 
function has previously been considered from a RMT-perspective in Refs. [15, [IB] and 
semiclassically in diagonal approximation [HI [HI [18] and beyond [19] . It is defined in 
terms of Z{6) given in (|8]) as 



A(.)=(JL 



- 



NJ 
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N 



' 2 77 S I ) 

1 N 



(24) 



n=0 



In particular, A(s) is the generating function for the variance (|a n | 2 ) of the coefficients 
of the characteristic polynomial. Note that |a n | 2 = \a^ 
function. In terms of sub-determinants we find 

(Kl 2 ) = E ( d rd^')u ■ 
r,r':[r[=|r'[=n 

We will show below that this reduces to the diagonal sum r = V for some specific 
ensembles. 



2 ensures that A(s) is a real 

(25) 



3.1.2. Parametric cross correlation for the spectral determinant Let U be a fixed 
unitary matrix and define U T := e lTPv U where r is a real parameter and P v is the 
projector onto the v-th basis state; the corresponding matrix is zero everywhere apart 
from one unit entry at the f-th diagonal position. Physically one may think of the 
parameter variation of a local magnetic field. Denoting the corresponding 

coefficients of the characteristic polynomial as 



r : |r|=n 



T:\T\=n, 



we consider the following parametric correlation function for the spectral determinant: 

1 

N 



B{t) 



d9 Z U(T) (9) Z v (9) 



u 



-ir/2 



E Mr)a n (0)*} 



u 



(27) 



n=0 



The above expression reduces the problem to the parametric correlations of the 
coefficients a n (r) which can be expressed as 



/ 



(a n (r)a n (0)*) l 



E 

r':|r'|=n 



\ 



(drd* r ,) + e iT ^ (d T d* r ,) 



r : |r|=n, 



r:|r|=7? 



(2f 



/ 



The two inner sums are here restricted to sets T of size |T| = n such that the marked 
v-th basis state is not in V for the first inner sum and the marked basis state is an 
element of T for the second inner sum. The outer sum over T' is only restricted by 



n. 



For the quantity B{r) we can now reduce the number of terms ton < N/2 by using 
the relation d r (r) = e l ^ +T ^ d*Jr) , Eq. (11). Note that the relation for the a n in Eq. rt5) 



would not be sufficient here, as the different components contributing to a n are exposed 
to different magnetic fields. 

As a side remark let us mention here that one can give an independent proof of the 



identity (11) in 2.2 using the reality of expressions of the type Zu{ T ) if one replaces P v 
in U T W = e lrPv by a unitary diagonal matrix with iV independent diagonal entries and 
expands the expression into a multinomial. 
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3.1.3. The spectral two-point correlation function and the form factor The spectral 
two-point correlation function is defined as 

R 2 (s) :=A 2 (J- d6p(6 + As/2)p(6-As/2)^ -1 (29) 

where A = j^Xfl)) = 2n/N is the mean spacing between eigenphases. 
Expanding the density of states in terms of traces and performing the integral over 
8, one obtains the standard expression 



, , 2 *^ / 27m \ 
^00= A? Z) 008 (a— )K n , (30) 



N ^ \ N 

n=l 



where 



Kn=±(\tTU*\X (31) 

is known as the form factor. The form factor played an important role in understanding 
universal and non-universal aspects of spectral statistics; here we give a new 
representation in terms of sub-determinants, that is, 

K n = ^ £ ( (m-1! m'-UlI j — Tj~^ — 

JM \ ±J - mAm-l 

m,m':|mr| = |mT|=n\ j=l •> i jj 

(32) 

This is an exact expression for the form factor for any ensemble of unitary matrices. 
We will show below that for some standard models, the double sum over multiplicities 
m and m' can be restricted further. 



3.2. Random-matrix theory 

Let us now consider unitary N x N matrices U which are distributed according to the 
Circular Unitary Ensemble (CUE) - in other words U has a uniform distribution with 
respect to the Haar-measure on the unitary group U(N). The spectral fluctuations of 
this ensemble are very well understood with explicit results for a large number of relevant 
measures. These known results have many implications for the statistical properties of 
the sub-determinants. 

One obtains, for instance, for the correlations of the coefficients a n of the 
characteristic polynomial [15] 

{a n a* n ,) cue = C' ( a n<v)cuE = ; (33) 

it is straight forward to extend this result to the correlations between sub-determinants. 
Indeed, any average over CUE is necessarily invariant with respect to conjugation, left 
multiplication, and right multiplication, that is, U i— > VUV' ,VU,UV with a unitary 
matrix V. Choosing V diagonal implies 

(drdr)cvE = (34) 
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and 



Note that 



{d r dr,) C xjE = d~rr'Cr ■ (35) 

22 c r = (K| 2 )cue = 1 (36) 
r : |r|=n 

where the sum extends over ( ) contributions. Moreover, if T and V have the same 
size, that is, |T| = |r'| = n, invariance of the ensemble average under conjugation with 
a permutation matrix implies 

c r = c r = c n c n = ( J . (37) 



Let us now consider the parametric correlation B(r) defined in (27). Note that it will 
not depend on the marked basis state, as the double sum over T and V in Eq. p8| only 
contains diagonal expressions after the CUE-average. Moreover among the (^) subsets 
T of a given size |T| = n > there are (^Zi) subsets which contain the marked basis 
state and all give the same contribution such that 

ne tT + N — n 

w 

and 



(a n (r)<) CUE = (3£ 



5(r) CU E = (iV + 1) cos(r/2) . (39) 
Let us finally look at the form factor; the CUE result is 

{n 
— if n < N 
N (40) 
1 if n > N. 

We may compare this to the CUE average of the form factor expressed in terms of sub- 



determinants (32 ). Invariance of the CUE ensemble with respect to group multiplication 



and unitary conjugation restricts the double sum over multiplicities m and m' in (32). 
For example, invariance with respect to multiplication with diagonal unitary matrices 
implies that only those pairs can survive, for which the corresponding product of sub- 
determinants II^i" 1 an d n^i" 1 ^Tj v ^ s ^ eacn basis state with the same multiplicity; 
here, the multiplicity of a basis state is the number of times a given index appears in 
any pseudo-orbit of the product Ylj d™^ ■ Note that this does not imply rrij = m'- as 
there may be many choices for the multiplicities rrij of the subsets Tj that result in the 
same multiplicities of a basis state. Comparing the resulting expression with the exact 



CUE result (40) one may obtain a large set of identities that have to be obeyed by the 



correlations among the sub-determinants. 
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3.3. Quantum graphs 

3.3.1. Star graphs - an introduction A quantum graph is a model for a quantum particle 
that is confined to a metric graph. To keep the discussion simple we will only discuss star 
graphs which consist of one central vertex and N peripheral vertices. Each peripheral 
vertex is connected to the center by a bond (or edge) of finite length < Lb < oo. By 
L = diag(Li, . . . , Ljv) we denote the diagonal matrix that contains the lengths on its 
diagonal. On a given bond b we denote by Xb G (0, L&) the distance from the central 
vertex. A scalar wave function on the graph is a collection of N complex (square- 
integrable) functions &(x) = (^i(xi), . . . ,i/)n{ x n))- The wave function is required to 
solve the free stationary Schrodinger equation on each bond, at given energy E = k 2 . 
This implies ipb(%b) — &b (e lkXb + e ~ tkx b+^L b \ w j iere a& j s amplitude of the outgoing 
wave from the central vertex and we have imposed Neumann boundary conditions at 
the peripheral vertices (at Xb = Lb). The matching conditions at the central vertex are 
given in terms of a unitary N x N scattering matrix S which relates the amplitudes et& of 
outgoing waves to the amplitudes abt 2lkLh of incoming waves by a& = Sbb , e 2lkLv ay ■ 
Equivalently 

a=£f(Jfe)a (41) 
for the quantum map 

U{k) = T{k)S where T{k) = e 2ikL . (42) 



The condition (41) is only satisfied at discrete values of the wave number which form 
the (wave number) spectrum of the graph. As a side remark let us also note that the 
above defined quantum map for a star graph also describes the quantum evolution on 
directed graphs with first-order (Dirac-type) wave operators and bond lengths 2Lb [T7] . 
A more general quantum graph requires a description in terms of a 2iV x 2iV matrix 



Spectra of quantum graphs and spectra of the associate unitary quantum maps 
U(k) have formed a paradigm of quantum chaos due to the conceptual simplicity of 
the models. In fact, both types of spectra are to a large extent equivalent [20J and we 
will focus the present discussion on the spectrum of the quantum map U(k). It can be 
considered as an ensemble of unitary matrices parametrised by k. The corresponding 
average will be denoted by 

1 ' K 



(F(U(k))) h = Km - / dkF(U(k)) . (43) 

Note that the wave number k enters the quantum map U(k) = T(k)S only through the 
diagonal factor T(k) = e 2lLk . 

The sets T C X in this model are one-to-one related to sub-graphs spanned by the 
corresponding bonds. The sub-determinants dr of U(k) can thus be written as 

d r = e ikLr d r (44) 

where L T = 2 Xlfeer is twice the metric length of the sub-graph connected to T 
and d r = det Sr is the corresponding sub-determinant of the scattering matrix S. 
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A generic choice of lengths L& implies that the lengths are rationally independent 
(incommensurate), which will be assumed in the following. Incommensurability implies 



that ( e ^£&=i™& L 



> does vanish except for nib 

k 



for all b 



N. 



3.3.2. Results for general star graphs It is straight forward to implement the averages 
for the spectral fluctuation measures introduced in Sec. |3.1[ Let us start with the 



variance of the coefficients of the characteristic polynomial, Eq. (25), which build up 
the auto-correlation function A(s). Due to the difference in the metric lengths of the 



corresponding sub-graphs only diagonal entries in the double sum of Eq. (25) survive 
the average, that is, 

(M 2 ) fe = E Kf • ( 45 ) 

r : |r|=n 

Note that the expression can not reduce further due to averaging. Contributions from 
different sets T contain orbits of different length, so non-diagonal contributions made 
up of products of orbits from different sub-graphs T do not survive the average; orbits 
and pseudo-orbits contained in dr cover the same sub-graph T, and thus have all the 
same lengths [21] • In full analogy, we find 



(a n (r)a n (0)*) k = \d T 



r:|r|=«, 
v<£T 



r : |r|=n 



(46) 



for the parametric correlations (28). In contrast to the CUE result this will generally 



depend on the marked v-th basis state. 

Furthermore, for the spectral two-point correlations, the form factor reduces to 

2^-1 



n 
N 



E (H-i)!(im'i-i: 



L^LZn m,m / 

L m r=L m 



n 



-dr 



y(-d* r ) m i 



mjirrij 



(47) 

where the C n is the set of all lengths that are a sum of n (not necessarily different) bond 
lengths of the graph. We have used the short-hand notation L mr = J2j f =i m j^ j r j - Note 
that equality of metric length L m r = L m >r implies equality of the topological length 
|mr| = |mT| while the opposite is not true. Eq. (47) expresses the form factor as a sum 



over all possible metric lengths with a fixed number n of bonds and a sum over pairs of 
completely reduced pseudo-orbits of topological length n of the same metric length. 



3.3.3. The two-star graph It is instructive to work out the simplest non-trivial case 
iV = 2 in more detail. In this case the only choices for V are the empty set, 1^ = {1}, 
r 2 = {2} and T 3 = {1,2} with lengths L Fl = 2L h Lp 2 = 2L 2 , and L F:j = 2{L X + L 2 ), 
see Fig. [T] The zeta function can be described in terms of the sub-determinant 
(in = 5ii which is just the reflection amplitude from the first bond and the determinant 
dr 3 = det S = alone; without loss of generality we set det 5 = 1. The remaining 
relevant sub-determinant is given by dr 2 = S 2 2 = due to (11) and det S = 1. 
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Figure 1. The 2-orbit graph consists of one vertex and two bonds labeled 1 and 2. 



Let us consider how an expansion of the zeta function in terms of periodic orbits and 
pseudo-orbits as discussed in Sec. 2.3 would look like. By expanding the product (17) 
and reordering the terms according to a cycle expansion [2], one obtains, for example, 



(u(k)(0) 



(t 



12 



tit 2 )e 

-3i9 



l-(h + t 2 )e- 1 

— (tll2 — ^1^12 + ^122 — t\ 2 t 2 )e 

Writing this out in terms of determinants yields instead 



-~2i6 



ClJ (h)(0) 



A2L\k-i9, 



,i2(L 1 +L 2 )k-2ie 



(4£ 



(49) 



1 d {l} - e l2h2k - M d\ l} + e 
= 2e l{Ll+L2)k ~ ie (cos(2ife(Li + L 2 ) - 0) - K (e i{Ll - L2)k d {l} 

The cancellation of contributions from longer pseudo-orbits \p\ > 2 appearing in 
the expansion, Eq. (48), becomes obvious when writing periodic orbits as completely 
reduced pseudo-orbits. For example, the contribution ti 22 from the orbit {122} is exactly 
canceled by ti 2 t 2 from the pseudo-orbit {12}{2} contributing just with opposite sign. 
By applying this cancellation mechanism recursively, i.e. reducing the orbits step by 
step, also the cancellation of contributions from longer pseudo-orbits can be understood. 
A similar cancellation argument is also used by the cycle expansion. This is however 
different for the contributions t\ 2 — t\t 2 in (48). In this case a reduction of the connected 
orbit leading to cancellation is not possible. 

The equivalence between pseudo-orbits on a subset Y and its complements can be 
made more explicit. The first and the last term in Eq. (49) resulting from pseudo-orbits 
of zero length and the length of the full graph, respectively, both have modulus of order 
1 and yield a real contribution to (u(k)(0) when the phase factor e l ( L i+ L ^) k - te [ s taken 
out. The same holds for the second and the third contributions to Eq. (49) from the 
orbits on the set T = {1} and T = {2}, respectively. Here, the identity (11 ) comes in to 
yield a real contribution (up to an overall pre-factor). 

For this simple example, we can calculate the spectral measures discussed in Sec. 



3.1 explicitly. For the auto-correlation function, Eq. (24), one obtains 
A(s) = 2cos(vrs) + 2|J ri | 2 • 



For the parametric correlation function, Eq. (27), we consider U(k; r) 



(50) 

diag(e ir , l)U(k). 
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One then obtains 

B(t) = 2(l + \d Tl \ 2 ) cos(r/2) . (51) 



Eventually, let us consider the form factor K n for a given n as presented in (47). It 
contains a sum over pairs of multiplicities m = (mi, 7712,7713) and m' = (m^, m' 2 , m 3 ). 
Both sums are restricted to have the same topological length |mr| = |mT| = n which 
implies two restrictions, namely ra\ + m 2 + 2m 3 = n = m[ + m' 2 + 2m' 3 . Furthermore 
only pairs of multiplicities contribute that have the same metric length L m r = L m T or 
L\(mi— Tti'i+mz— m' 3 )-\-L2(m2— m' 2 +m 3 — m' 3 ) = 0. The latter implies m\ +m 3 = m^+m^ 
and m 2 + m 3 = ml 2 + m' 3 . Only three of these four restrictions on pairs of orbits are 
independent. The form factor can then be written as 

Kn = ~0 Srn 2 +m 3 ,m' 2 +m' 3 (n ~ m 3 - l)\(n - m' 3 - 1)! X 

0<m3 <ra/2 0<m2 <7i — 2tti3 
0<m f 3 <n/2 0<m^<n-2m^ 



n - 2m 3 \ fn - 2m! 3 \ (-l) m 3+m 3 |J ri |2(™-m3-m 3 ) 
m 2 J\ m'2 J m%\m 3 \(n — 2m 3 )\(n — 2m' 3 )\ 
Writing the Kronecker as 

1 r 2n 

S , / , / = / do p ia (™2- m 2+™3- m 3) 

makes it possible to sum over m 2 and m' 2 independently. With 

n-2m3 



(52) 



J- e Mrn 2 +ms)^n 2m^ =(2cos(a/2)) 



n— 27713 gicm/2 



and 



m,2=0 



n— 2m 3 



y- 3 e - ia(m2+m3 )^ 2m(A = (2cos(a/2))"- 2m 3 e - i " i! 2 

m' 2 =0 V 27 



we obtain 



n 2 f w sr^ fn- m 3 \ fn- m' 3 \ 



n - m 3 \ /n - m\\ (-l) m3+m sU r |2(n-m 3 -m 3 ) 

-if " x 



0<m 3 <n/2 
0<m 3 <n/2 



(77 — m 3 ) (n — m' 3 ) 



a \ 2(n-m 3 -m 3 ) 

2 cos - ) . (53) 



2 

The sums with respect to m 3 and m 3 can be performed by using 



0<m 3 <n/2 



/o (tj — 777.3) \ m 3 J 2 n n LV / 

(l - Vl -4:X 2 Y . (54) 
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This yields for K n 



Kr, 



1 

47T 



2?r 



da 



a 



cos 



d 



2n 



+ \ COS z 



a 



dr-. 



1 + 







/ 2 " 






loos- 




— A / COS z — 

V 2 




■- 



2/i 



(55) 



The constant term at the end describes the behavior for n 3> 1, the other two 
contributions describe oscillations around the asymptotic value K n = 1. By taking 
into account that the arguments of the square roots above are negative, we can rewrite 
the expression as 



K n = l + 



2vr 



2tt 



da cos 



2n arccos 



a: 



cos 



(56) 



The expression in Eq. (56), which we obtained from periodic-orbit expansions, coincides 
for all dr 1 with the result obtained in [23j starting from the eigenvalues of the quantum 
scattering map. This is the first derivation of the result (56) from periodic-orbit 
expressions for general dr l - In [23] 



was only done for 



d\ 



1/V2. 



a derivation based on periodic-orbit expressions 
we also take into account contributions beyond the 



Note that in contrast to 
diagonal approximation. Due to the factor (— l) m 3+«^ appearing in the Eqs. (52 53) 
the ones with m 3 — m' 3 odd contribute with negative signs leading together with the ones 
with mz — m' 3 even to a form factor smaller than expected in diagonal approximation; 
as expected it tends to K = 1 for n — > oo. 



4. Conclusions 



The goal of this article is two-fold: first of all, we advocate considering sub-determinant 
expansions for spectral functions and statistical measures such the density of states or 
various correlation functions. This makes it possible to separate out contributions which 
vanish after averaging and those whose non-diagonal contributions survive averaging. 
Secondly, we introduced a sub-determinant identity due to unitarity which makes it 
possible to give detailed relations between short and long orbits on a graph. In 
particular, this identity implies that contributions to the characteristic polynomial 
originating from irreducible pseudo-orbits of a given sub-graph have the same weight as 
the irreducible pseudo-orbits of the complement of that sub-graph and are additionally 
linked through a common phase factor. The identity leads to simplified expressions 
for the characteristic polynomial, the Newton identities and the spectral density. 
Furthermore we study the effect of this identity on spectral correlation functions such 
as the auto-correlation function of the characteristic polynomial, the parametric cross 
correlation function and the spectral form factor. 

We derive explicit expressions using sub-determinant expansions for a simple model, 
star graphs consisting of N bonds connected by a single vertex. We then work out 
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in more detail the simplest case N = 2. The new identity is essential to obtained 
the behaviour of correlation functions for small energy differences or large times. It 
captures additional correlations between orbits of different length and needs to be taken 
into account when singling out correlated orbits which survive averaging. 

In this context several potentially interesting extensions arise: taking the 
semiclassical limit on both sides of Eq. (22), the two expressions are semiclassically 
not obviously identical. The left-hand side leads to the Gutzwiller trace formula which 
contains orbits of arbitrary length while the right-hand side contains pseudo-orbits of 
finite length (and their repetitions). For short orbits n <^ N, one may argue that the 
two expressions have a semiclassically small difference, for longer orbits this is far less 
obvious. 

A second point concerns the exponential proliferation of the number of orbits in the 
standard trace formulae. It is tamed to a certain degree when using sub-determinants 
by the fact that different contributions contribute with different signs. Thus the sub- 
determinant expressions contain large fluctuations. Understanding overall cancellations 
is an interesting task. For example, the form factor for the two-star graph for large n 
contains positive and negative contributions which on their own grow as n — > oo while 
their difference remains 0(1) as can be checked from the expressions for K n given above. 

The analysis of spectral correlations focused here on the general unitary case. It 
would be interesting to include the effect of self crossings of orbits that allow for partners 
traversing parts of the diagram in different directions. This would capture effects arising 
due to time reversal symmetry. 

For graphs the supersymmetry technique gives an alternative approach to obtain 
universal results [2S]- With supersymmetry one may derive universality under 
sufficiently nice conditions, however, rigorous proofs are still not available. The main 
obstacle for the supersymmetric approach seems to be repetitions which are difficult to 
incorporate correctly [26] . The proposed approach may help us understanding the effect 
repetitions on spectral correlation functions. 
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